Calculation of bulk viscosity of QCD at zero temperature and finite chemical potential 
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In this letter, based on Kubo's formula and the QCD low energy theorem, we propose a direct 
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Recently transport coefficients of hot/dense quark 
matter attract lots of attentions. It is observed that the 
quark gluon plasma (QGP) created at RHIC (Relativis- 
tic Heavy Ion Collider) behaves like a nearly perfect fluid 
[l|-y] with the ratio of the shear viscosity 77 over entropy 
density s approaching the lower bound of I/Att [5| near 
the phase transition temperature Tc- Such a behavior 
indicates that the created QGP is not a weakly coupled 
gas which is expected from perturbative quantum chro- 
modynamics (pQCD), but a strongly coupled liquid. On 
the other hand, lattice QCD calculation shows that the 
ratio of the bulk viscosity C, over entropy density s has 
a sharp peak near the critical temperature [y, lZ| • Such 
a sharp peak behavior of C, has also been observed in 
many model calculations including linear sigma model 
[8|, real scalar mod el [QJ , massless pion gas model [lOJ, 
quasiparticle model [ll[, Z(2 ) an d 0(4) models I4I, and 
Dyson-Schwinger equations [Ij]. The bulk viscosity C, 
characterizes the response of the system to the confor- 
mal transformation and vanishes when the system has a 
conformal equation of state, therefore the sharp peak of 
the bulk viscosity would strongly affect the physics of the 
QCD matter near Tc and is very important for the study 
of QCD phase transition. 

Because the phase structure of QCD at finite density 
is much richer than that at finite temperature (see, for 
example, Ref. [ij]), one may naturally ask whether the 
sharp peak behavior of C, exists at zero temperature and 
finite /i. As far as the present authors know, up to now 
there are only limited number of papers trying to answer 
the above question. For example, in Ref. 15| the authors 
study the viscosity at finite m with Nambu-Jona-Lasinio 
(NJL) model and in Ref. [16| the authors study the vis- 
cosity of strange quark matter at finite /x with quasi par- 
ticle model. Just as will be shown below, elucidating 
the bulk viscosity of QCD matter requires a knowledge 



of the equation of state (EOS) of strongly interacting 
matter, and the plausibility of the results of any given 
study rest on the EOS employed, which is difficult to 
judge a priori (Presently, it is not feasible to study the 
EOS at finite /i rigorously in lattice simulation of QCD, 
due to the problems arising from the complex valued Eu- 
clidean fermion determinant at nonzero chemical poten- 
tial). Consequently, the exploration of alternative equa- 
tions of state and the identification of qualitatively con- 
sistent results is important. Here we would like to stress 
that the absence of confinement is a defect common to 
most existing model studies for the bulk viscosity of QCD 
matter [sMlalla. Il6l|. A covariant, continuum approach 
incorporating both confinement and chiral symmetry in 
that direction is necessary. 

It is well known that the Dyson-Schwinger equa- 
tions (DSEs) provide a nonperturbative, Poincare invari- 
ant continuum framework for analyzing QCD, i.e., for 
the study of confinement and dynamical chiral symme- 
try breaking (DCSB), and hadron observables |17h21| . 
Herein we will employ a dynamical, confining Dyson- 
Schwinger equation model of QCD [2^ to calculate the 
bulk viscosity of QCD at finite chemical potential and 
zero temperature. 

In order to make this letter self-contained, let us first 
recall the general approach of calculating bulk viscosity 
in existing literature. According to Kubo's formula the 
bulk viscosity Q can be expressed as [6| 
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where the retarded Green's function G^ is defined as: 



G^{x) 



diod^k 



G^(cj,fc)e~*('^*~*^-^) 



(2^)4 



(2) 



with 7^1^ the stress tensor of the system. In the small fre- 
quency region one can assume the spectral density p{uj,p) 
which is defined as 
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where the parameter wq identifies the scale at which the 
perturbation theory becomes valid [6| . With such ansatz 
the bulk viscosity can be expressed as 
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For QCD, the trace of stress tensor reads 
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where g is the strong coupling constant, Tp and To are 
the contribution of quark fields and of gluon field, re- 
spectively, and (3{g) is the QCD /3-function which deter- 
mines the running behavior oi g. In Eq. ([5]) q are quark 
fields with two fiavors (in this letter we will limit our- 
selves in two flavor case and set the current quark mass 
TTiu = rud = m). From the QCD low-energy theorems at 
finite temperature T and /i [23j] , one can find 

^^^^Tp' rf) {o)t = j d^x{T,{TG{x)dmil) 

where d is the canonical dimension of the operator O. 
Using the above equation, one has 
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From the above two relations one obtains 
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where £ is the energy density and V is the pressure den- 
sity of QCD. Here, because the current quark mass m of 
u and d quark is very small, in deriving Eq. (jlOp we have 
neglected the term proportional to m? . 

Now we take the T — >■ limit of the above expression. 
In this limit the terms involving Td/dT vanish. There- 
fore one would obtain the following expression of C, at 
zero T and finite p: 

%Uo = -^pp^ip) + l(^V[p) + 2mp^M. _ Am{qqUl) 

where Pq{p) is the quark number density which can be 
written as (see, for example, Ref. [24]) 
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and {qq) is the quark condensate at zero T and finite p 
which can be defined as 

^tr{G,[^](p)}. (13) 



In Eqs. ([T2l) and (fT3|) Gq[p]{p) is the dressed quark prop- 
agator at finite /i, Nc and Nf denote the number of colors 
and of fiavors, respectively, and the trace operation tr is 
over Dirac indices (we will not write the renormalization 
constants explicitly because one would find in the final 
result the renormalization constants cancel each other). 

Eq. (ITT|) indicates that the bulk viscosity of QCD at 
zero T and finite p is determined by Pq{p), the pressure 
density Vip) and the quark condensate at finite p. From 
Eqs. ([121 and (TTSl) it can be seen that the quark number 
density and quark condensate at finite p is totally deter- 
mined by the dressed quark propagator at finite p. As 
will be shown below, the pressure density is also totally 
determined by the dressed quark propagator at finite p. 
This means that the bulk viscosity at finite p is totally 
determined by the dressed quark propagator at finite p. 
This is quite different from the case of finite T in which 
the term involving the speed of sound contributes to the 
peak behavior of (/s around Tc [6|. 

From Ref. [2J] one can obtain the following: 
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Based on the above equation it can be seen that the pres- 
sure density Vip) is the sum of two terms: the first term 
2 ) {Tf)t V{p = 0) (the pressure density of vacuum) is only a 
/i-independent term; the second term, which is totally 
determined by Gq[p]{p), contains all the nontrivial p- 
dependence. Just as was shown in Ref. ^24], the con- 
, stant term V{p = 0) can be calculated through Cornwall- 
^ ' Jackiw-Tomboulis effective action [25[. However, in the 



present paper we are only interested in the ^-dependence 
of C, so we ignore this term. 

According to Eqs. (8-10) and Eq. (fT4|) . once the 
dressed quark propagator at finite /^ is known, one can 
obtain the bulk viscosity of QCD at finite ^. However, at 
present it is very difficult to calculate Gq [/i] [p) from first 
principles of QCD and one has to resort to various QCD 
models. Over the past few years, considerable progress 
has been made in the framework of the rainbow-ladder 



approximation of the DSEs approach 17H2l| , which pro- 
vides a successful description of various nonperturbative 
aspects of strong interaction physics. In the rest of this 
paper we will use DSEs approach to calculate C- 

To obtain the dressed quark propagator at finite /i let 
us first recall the following general result given in Refs. 
|26l - l3q | : If one adopts the rainbow approximation of the 
DSEs (rainbow approximation means one uses the bare 
vertex to approximate the full quark-gluon vertex in the 
quark DSE) and ignores the /x dependence of the dressed 
gluon propagator (one often makes this approximation 
in the study of finite density QCD because it is usually 
believed that the effect of chemical potential on the gluon 
propagator arising from quark-loop insertions is small in 
comparison with that on the quark propagator; this is 
a commonly used approximation in studying the dressed 
quark propagator at finite /x [18J, I26l435l |). then one can 
obtain the dressed quark propagator at finite jjl from the 
one at /i = by the follo wing substitution (more details 
can be found in Refs. [26l43d |): 
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where p = {p,p4 + ifi). The validity of this approximation 
has been discussed in detail in Refs. 28143d. l36l | and has 
beer i ap plied to the study of properties of quark stars 
371 138[. Therefore, once the dressed quark propagator 



at /i = is known, the dressed quark propagator at finite 
/i and r = can be obtained by means of Eq. (TTSl) . 
Now one needs to specify the form of the dressed quark 

I 



propagator at ^ = T = 0. In Ref. [22|, guided by the 
solution of the coupled set of DSEs for the ghost, gluon 
and quark propagator in the Landau gauge, the follow- 
ing meromorphic form of the dressed quark propagator 
is proposed: 






ip + aj + ibj 



i^+ aj — ibj 



(16) 



The propagator of this form has np pairs of complex 
conjugate poles located at aj ± ibj. When some bj is 
set to zero, the pair of complex conjugate poles degener- 
ates to a real pole. The residues tj are real (note that a 
similar meromorphic form of tjie quark propagator was 
previously proposed in Ref. ,3^], in which the residues 
in the two additive terms are complex conjugate of each 
other) . The values of parameters in Eq. ^TE\i are listed 
in Table I 1221 where 2CC means two pairs of complex 
conjugate poles, IRICC means one real pole and one 
pair of complex conjugate poles and 3R means three real 
poles. The parameters of these meromorphic forms in 
Table I are fitted by requiring that these meromorphic 
forms of propagator reproduce some physical observables 
(the two-quark condensate and pion decay constant) and 
lattice data (the zero momentum values of the mass func- 
tion and wave-function renormalization, Mq and Zq , and 
an approximate width of the region of large dynamical 
mass generation, ujl) (for more details, please see [22|). 
Here we should stress that the real dressed quark prop- 
agator has much more complicated analytic structures 
than that given in Eq. (J16p . However, despite its sim- 
plicity, this kind of parametrization of the dressed quark 
propagator has several appealing features: it can simulta- 
neously describe DCSB and quark confinement which are 
very important for low energy QCD; at large momenta 
it reduces to the free fermion propagator, as required by 
asymptotic freedom. Phenomenologically, these param- 
eterizations of the dressed quark propagatoi^ can give a 
quite good description of low energy QCD [22j |. 



Tabic I. The parameters for the quark propagator. These parameters are taken directly from Table II of Ref. |22ll . 



Parameterization 


n 


ai (GcV) 


bi (GeV) 


r2 


02 (GeV) 


fc2 (GeV) 


ra 


aa (GcV) 


2CC 


0.360 


0.351 


0.08 


0.140 


-0.899 


0.463 


- 


- 


IRICC 


0.354 


0.377 


- 


0.146 


-0.91 


0.45 


- 


- 


3R 


0.365 


0.341 


- 


1.2 


-1.31 


- 


-1.06 


-1.40 



With such propagator the quark number density and 
the M-de pen dent term of P{ii) in Eq. (fT4)) can be calcu- 
lated as [24 1 
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Froni Eq. ([TT]) it can be seen that when n is below a 
critical value fic = ^^{Wjl} (a'c = 351 MeV, 377 MeV 
and 341 MeV for the 2CC, IRICC and 3R parametriza- 
tion, respectively), the quark number density vanishes 
identically. Namely, /i = /ic is a singularity which sepa- 
rates two regions with different quark number densities. 
This result agrees qualitatively with the general conclu- 
sion of Ref. [40|. In that reference, based on a universal 
argument, it is pointed out that the existence of some sin- 
gularity at the point n = n^ and T = is a robust and 
model-independent prediction. The numerical value of 
the critical chemical potential in pure QCD (i.e., with 
electromagnetic interaction being switched off) is esti- 
mated to be 307 MeV. The value of Hc obtained here is 
of the same order of magnitude as the estimates in Ref. 

[li. 

The quark condensate at finite // with dressed quark 
propagator given in Eq. ([TC|) has been obtained in Ref. 
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where {qq)o is the quark condensate ed, fi — and fj 
reads 
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From Eq. ([^0]) it can also be seen that when ^ < ^c 
quark condensate equals its vacuum value. 

Substituting Eqs. (HZl), ^ and 1^ into Eq. ^ 
one can calculate the bulk viscosity at zero T and finite 
fi. The numerical result is shown in Fig. 1. From Fig. 
1 one can find the bulk viscosity is kept unchanged from 
the value at fi = when fj, is smaller than fic^ and fic is 
the point when quarks get populated (like the liquid-gas 
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FIG. 1. Bulk viscosity at finite fi 



point, if there were a self-bound nuclear or constituent- 
quark groundstate). This is understandable because as 
was analyzed in Ref. |24l |. when fi < n^, the ground 
state at finite /i remains the vacuum and the physical 
quantities should be kept unchanged. When fi > fj,c, 
it can be seen from Fig. 1 that the bulk viscosity de- 
creases monotonously as n increases and no sharp peak 
is observed at finite fi. This behavior of bulk viscosity at 
finite fi does not show a first-order chiral phase transi- 
tion (or a rapid cross-over) , which is quite different from 
that at finite temperature. The result of such a behav- 
ior is qualitatively consistent with the result obtained in 
Ref. [15J. Here, we would like to point out that the 
NJL model [4^ |43| is the most widely adopted nonper- 
turbative QCD model in the study of finite density QCD 
at present. The NJL model can phenomenologically in- 
corporate DCSB but does not have quark confinement, 
whereas the DSEs approach can simultaneously incorpo- 
rate DCSB and quark confinement. Over the past few 
years considerable progress has been made in the frame- 
work of the rainbow-ladder approximation of the DSEs 
approach at finite temperature and finite density (see, 
e.g., the review article Ref. 18]). Due to this we natu- 
rally expect that the DSEs approach can be applied to 
the study of bulk viscosity of QCD. This is our main mo- 
tivation to use the DSEs approach to calculate the bulk 
viscosity of QCD at finite density in this work. Finally, it 
should be noted that that in obtaining the result of bulk 
viscosity in this paper, we have made use of Eq. (1151) . 
However, 



as was pointed out in Refs. [26l - l30j . Eq. p5)) 
is obtained under the approximation that one neglects 
the fj, dependence of the dressed gluon propagator. This 
approximation is valid only in small fj, region. Therefore, 
if one applies Eq. ([T5|) to the case of large fi, one should 
be cautious. 

To summarize, in this letter, based on the Kubo's for- 
mula and the QCD low energy theorem, we have derived 
a direct formula for calculating the bulk viscosity of QCD 
at finite chemical potential /i and zero temperature. Ac- 
cording to this formula, the bulk viscosity at finite fi is 
totally determined by the dressed quark propagator at fi- 
nite fj,. Then, under the rainbow approximation of DSEs 
and ignoring the /.t-dependence of the dressed gluon prop- 



agator, we use the meromorphic quark propagator pro- 
posed in Ref. 22| to calculate the bulk viscosity at finite 
/i. It is found that when ji is below a critical value ^ci 
the bulk viscosity equals its vacuum value. When ^ > ^^j 
the bulk viscosity decreases monotonously as /i increases, 
and no sharp peak is observed at finite /i. This means 
that in our model calculation in the framework of DESs 
approach the bulk viscosity at finite ^ does not show a 
first-order chiral phase transition (or a rapid cross-over) , 
which is quite different from that at finite temperature. 
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